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Abstract 

We consider one-dimensional classical time-dependent Hamiltonian systems with quasi-periodic 
orbits. It is well-known that such systems possess an adiabatic invariant which coincides with the 
action variable of the Hamiltonian formalism. We present a new proof of the adiabatic invariance 
of this quantity and illustrate our arguments by means of explicit calculations for the harmonic 
oscillator. 

The new proof makes essential use of the Hamiltonian formalism. The key step is the introduc- 
tion of a slowly-varying quantity closely related to the action variable. This new quantity arises 
naturally within the Hamiltonian framework as follows: a canonical transformation is first per- 
formed to convert the system to action-angle coordinates; then the new quantity is constructed as 
an action integral (effectively a new action variable) using the new coordinates. The integration 
required for this construction provides, in a natural way, the averaging procedure introduced in 
other proofs, though here it is an average in phase space rather than over time. 

PACS: 45.05,+x, 45.20. Jj. 



1 Introduction 

At the Solvay Conference in 1911, Einstein answered a question raised by Lorentz with the statement 'If 
the length of a pendulum is changed infinitely slowly, its energy remains equal to hv if it was originally 
hv\ 1 The expression 'adiabatic invariant' came to be used in this context, and this has stuck over the 
years, even though the change of the system has to be slow, rather than adiabatic in the thermodynamic 
sense 2 and the quantity in question is not invariant — and indeed may not even change very slowly, as 
will be discussed further. 

Adiabatic invariance is readily understandable at a superficial level: students studying classical 
mechanics have no difficulty with the idea that, for a simple pendulum, the quantity E/u> changes 
very slowly in some sense when the frequency to is slowly varied. However, the details are often found 
more elusive, and with good reason: the concept is fundamentally subtle and the proof of invariance 
is genuinely difficult. The treatment given in most texts relies on the theory of action- angle variables, 
which already involves a significant level of sophistication; and in addition some form of averaging is 



1 This was apparently first noticed by Rayleigh in 1902 for the classical pendulum; however, Lorentz's question related 
to a pendulum oscillating in a particular quantum state. 

2 An adiabatic change in thermodynamics is one that happens with no gain or loss of heat. 
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usually required. Moreover, many texts are careless of details, omit to state exactly what is being 
done and do not reveal the motivation for doing what is done. There is a danger too, always present 
when canonical transformations are used, of losing track of the relation between the dependent and 
independent variables. 

We confine our attention to the case of a one-dimensional system governed by a Hamiltonian 

H(p,q,\(t)), (1) 

where A is slowly varying in a sense that we will define precisely in sectional The quantity / defined 
by 

I-^fpdq (2) 

is an adiabatic invariant of this system. The integral is to be taken round the curve (assumed to be 
closed) in phase space given by H(p, q,X) = E at constant A. It represents the area inside the curve. 
It is important to realize that this integral is a purely geometrical construction; it takes no account of 
the time evolution of the system. 

It is always a good idea to have in mind the example of the time- varying simple harmonic oscillator 
with Hamiltonian H(p, q,u>) given by 

H(p,q,oj(t)) = ^(p 2 +u J 2 (t)q 2 ). (3) 

For this system, we may easily calculate the integral in © to obtain I = E/u>, with the implication 
that to a good approximation the energy will change proportionally to the frequency. 

For many simple systems, it is not even necessary to calculate the integral. For example, in the case 
H = (p 2 + Ag 4 )/2, we have 

I = — (f y / 2E - Xq i dq= -(2E)i\-i [ y/l - x A dx (4) 
2vr / 7T y_ 1 

from which we see that J is a numerical multiple of (E 3 / A) * and that in this case to good approximation 
the energy varies as A a for slowly varying changes in A. 

Of course, I is not in general exactly invariant. 3 It is necessary therefore to give a precise definition 
of adiabatic invariance. Clearly, any quantity such as I (which is essentially an average over a cycle) 
will vary slowly, since the system varies slowly. In order for I to be adiabatic, it must vary more slowly, 
in some sense than A varies and this means that a proof of adiabatic invariance must keep careful track 
of any error terms. 

In the next section, we illustrate some of the features of I by examining a particular adiabatic 
variation of the harmonic oscillator. In the third section we discuss strategies for differentiating I in 
order to determine its rate of change. In the fourth section we revisit the harmonic oscillator problem, 
in this case we consider a general adiabatic variation and present a new proof of the adiabatic invariance 
of I. In the fifth section, we show how this proof can be applied to more general systems. Finally in 
section 6 we summarize our discussion and present our conclusions. 



2 The time dependence of /: harmonic oscillator case 

In some cases, it is possible to investigate the time dependence of I by solving Hamilton's equations 
and evaluating the area integral (0) explicitly. It is instructive to look at the simplest case of the simple 

3 Interestingly, for the harmonic oscillator, there is an exact invariant (called the Ermakov-Lewis invariant — see [Q 
for a useful discussion) which can be written in the form 

l((q/pf + (p q -qpf) 

where p is any solution of p + u) 2 p = p~ 3 . 



C.G. Wells and S.T.C. Siklos 



3 



harmonic oscillator © with 

lu 2 = 1 + et. (5) 

This could be regarded as the approximation, by Taylor series, to a more general slowly- varying u>(t), 
but for present purposes that would sacrifice clarity to little purpose: it would entail keeping track 
of an additional set of error terms and the calculations are messy enough in the special case. For 
similar reasons we will restrict attention to e > 0. The final expressions would look better if we had set 
uj 2 = u> 2 (l+et), but the necessary factors can easily be inserted at the end on dimensional considerations. 

We will solve Hamilton's equations explicitly for this system in terms of Airy functions, then calculate 
the value of I(t) using standard asymptotic expansions for the Airy functions. 

Hamilton's equations are q = p, p = —iu 2 q, leading to 

?+(l + et)g = 0. (6) 

Setting z = — e~ 2 / 3 (l + et) reduces this to the Airy equation q" — zq = with solutions Ai(z) and Bi(z) 
(see for details). Making use of the fact that the Wronskian of Ai(z) and Bi(z) is tt~ 1 , we can write 
the solution that satisfies q = and q = 1 at t = in the form 

Tr" 1 ^*) = Ai ( - e-^cj 2 (f)) Bi' ( - e -t) - Bi ( - e-%uj 2 {t)) Ai' ( - e"* ). (7) 

For Kl and et ^ 1, we can apply the standard asymptotic approximations to obtain 

q(t) =LU-i cos 9 + ^eoj~i sin 8 + 0{e 2 ); p(t) = -w% sin 9 + 0(e 2 ), (8) 

where 

Note that 9, which arises naturally from the asymptotic expressions, is well behaved in the limit e — > 0. 
We can now determine the time dependence of /: 

I = H/oj = ^{p 2 /u + uq 2 ) = - + iew~ 3 sin(9cos(9 + 0(e 2 ). (10) 

From this result, we can make two important observations: 

• The time rate of change of I(t) is not smaller than that of the Hamiltonian; both are first order 
in the small parameter e; 

• I(t) does not grow at a rate proportional to e; rather, it oscillates and the time-average does not 
contain a term proportional to e. 



3 The time dependence of /: general case 

In order to evaluate the integral (J2J), we can in principle solve the equation H(p, q,t) — E at any fixed 
time t, to obtain a solution of the form p = P{E,q 1 t). Of course, in practice the appropriate value 
of E at time t can be obtained explicitly only by solving Hamilton's equations for the system with 
appropriate initial conditions, obtaining a trajectory in phase space of the from p — P(t), q = Q(t) and 
setting E(t) = H(P(t),Q(t), t). 

At this point, it is instructive to analyse a calculation that is used the basis of the proof of adiabatic 
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invariance in a number of standard text books. We start by differentiating under the integral sign 
^fp(E,q,t) dq = 






1 dH 




dP 






dp 


-\ 

q,t 


' ~dt 





dq 



dP 




dP 




~dt 


-\ 

q.E 







dq = 0. 



(11) 
(12) 
(13) 
(14) 



This seems to be a good result, but as it stands here it is wrong. The step that is incorrect depends on 
how the integral is to be interpreted and three different interpretations can be found. It should be said 
that the proofs of adiabatic invariance normally include some sort of averaging which can obfuscate the 
deficiencies of the above argument. 

If the integral is interpreted as being taken round the closed curve H (p, q,t) — E (which is how 
such integrals will be defined in this paper), then equation l|12|) is wrong: the identity dE/dt = dH/dt 
follows from Hamilton's equations and Hamilton's equations for the trajectory of the system whose rate 
of change of energy is dE/dt do not hold on the closed curve. Instead, dE/dt = on this curve, which 
leads to a rather unhelpful integral. 

If instead the path of integration is along the trajectory of the system between points at which q = 
(say) , then the curve is not closed and attention would have to be paid to the variation of the endpoints 
in equation (|11J) : to put it another way, the differentiation cannot simply be taken under the integral 
sign. As a complicating factor in this case, the integral we are differentiating is not equal to J, and 
extra error terms need to be included to account for this discrepancy. 

Finally, if the path of integration is the closed curve E — H(p, q, t) at a fixed time t and then each 
point on the curve moves along the Hamiltonian flow to obtain a closed curve at later times, then the 
rate of change of the integral is indeed identically zero (by Liouville's theorem) but the later curve is no 
longer of the form E = H(j>, q, t) (for example, it does not remain elliptical in the case of the harmonic 
oscillator). 

The remaining equations are standard results: Ijl3(l comes from differentiating E = H(P(E, q, t), q, t) 
partially with respect to E and (fT4*jl comes from differentiating this same equation partially with respect 
to t. 

Although we do not use the above calculation explicitly in the proof of adiabatic invariance of /, it 
does provide a useful idea that underlies our method of proof. It is relatively easy to calculate / in any 
given case, but it is not at all easy to determine its time dependence. Conversely, there are quantities 
defined by integrals closely related to / which are not easy to calculate, but the time-dependence of 
which can be determined. The underlying idea of our proof is to define a quantity J that can be shown 
to be close to / and whose time dependence can be determined. 



4 Adiabatic invariance of /: the harmonic oscillator case 

The proof of the adiabatic invariance of / for a general Hamiltonian presented in the next section is 
quite difficult to understand at first because of the number of changes of variable and the different 
independent variables used. This is a problem common to all such proofs. Instead of starting with the 
full proof, it makes sense to see how it works in the more restricted case of the harmonic oscillator: 

H(p,q,T) = ±(p 2 +q 2 /T 2 ). (15) 

We have used r(t) as the slowly varying parameter rather than u because it has the dimensions of time. 
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As mentioned above, the plan is to construct a quantity J which remains close to / in value as time 
evolves, but whose time variation can be determined more easily. 

The essence of the proof is to write the system in action-angle coordinates (I, 4>) with the dynamics 
governed by the Hamiltonian K(I,4>,t). Since K is not independent of time we construct an action 
variable associated with this new Hamiltonian system. This new action variable J, is an average of I 
over the cycle defined by $S <fi ^ 2n. If the system were independent of time, then / and J would 
coincide. However when r is varying they may be shown to differ by an amount which is 0(f). The 
time evolution of J is much slower than that of / and may be determined from the fact that K may be 
written as a function of J and r alone (this follows from the fact that J is the action variable derived 
from the Hamiltonian K). In this way we are able to deduce that J = 0(f, f 2 ) and the adiabatic 
invariance of J and hence / may then be established. 

We start by defining a new pair of coordinates (I, (j>) by 

I=-(p 2 T + q 2 /T) = HT, (f> = sin' 1 —= q = sin' 1 - 9 . (16) 

I is the adiabatic invariant defined by the integral J5J. It may easily be verified that if p and q satisfy 
Hamilton's equations with respect to the Hamiltonian i|15[l . then / and <j) satisfy Hamilton's equations 
with respect to the Hamiltonian 

K(I,<f>,r) =I{l-i sin (p cos <j>)/r. (17) 
The coordinates (I, 4>) are therefore canonical 4 . 



Now we define the quantity J by 

1 



2tt 



J=-J q m (18) 

where the path of integration is the closed loop in the I,4> plane determined by K(I, 0, t) = constant 
and t = constant. Note that this is precisely analogous to the definition @ of /. Just as / can be 
written in terms of p, q and r, or equivalently in terms of H and r, we can write J in terms of /, <f> and 
t, or equivalently in terms of K and t, using equation (|17JI : 

J= K 1 [ 2 * d0 = ^^__^ > 

2-7T J 1 — f sin cos q> 

The value of the above integral can be most conveniently obtained using an integral round a unit circle 
in the complex plane. We can now relate J to / by again using equation (|17l) : 



J = 1(1 - f sin0cos0)v/l + f 2 /4. (20) 

Thus 

\J-I\=0(f) (21) 

which is the key result that ties / to J. 

Next we need to determine the evolution of J. Differentiating equation l|19fl gives 
dJ d / > nrrr\ / t^ttcLK H 



dt dt 

= T 



+ = 7Vl + f 2 /4— + K- (tVI + t-74) (22) 

v/TT7 ^(^L +a 1( tVTT7571 ) (23) 

v/l + r 2 /4^ (7(l-fsin0cos0)/r)+^^ (ry/l + f 2 /4) (24) 



|-) ((1 - f sin 4> cos 4>) y/l + i 2 /A) . (25) 



These coordinates may be derived from the type 2 generating function 

F(I,q,r) = J 9 pdq = J 9 ^ 21 /t - q 2 /r 2 dq = I sin' 1 (q/V^) + I(q/V2Ir) \/ 1 - q 2 /2Ii 
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Thus 

1 dJ — t shi (j> cos 4> tt 
J dt 1 - f sin 4> cos 4> 4 + f 2 ' 

This is the result we want. It shows that, if f and f are small, then J(t) will not vary much from 
its original value and nor, in view of (|2ip. will I(t). The details of this final step are deferred to the 
next section, once it has been determined exactly what is required from a formal definition of adiabatic 
invariance. 



5 Adiabatic invariance of /: the general case 

First we need to adopt a definition of adiabatic invariance. The most useful one for our purposes is as 
follows: 

Let T be an arbitrary fixed time, for all e > consider variations r(t) = t(c,<) such that r = 0(e) 
and t = 0(e 2 ) (e.g., r(t) = h(et) for some function h). We say that a quantity I(t) is an adiabatic 
invariant of the dynamical system if for all such variations of r we have \I(t) — I(0)\ = 0(e) for all 
O^K T/e, 5 

Now we follow the method laid out in the previous section for the special case of the harmonic 
oscillator. 

From the Hamiltonian H(p,q,T(t)) we construct new variables / and <j> as follows. First we define 

I ^hj >pdq ^hj >p{Hiq ' T)Aq (27) 

where the integral is taken over the curve (assumed to be closed) H (p, q, r) = constant at constant r. As 
before, / will emerge as a function of H and r. The generating function for the canonical transformation 
is given by 

F(I,q,T) = ppdq', (28) 
Jo 

where p in the integral is written in terms of /, q' (a dummy variable) and r using P7)l. Then the angle 
variable (f> and the new Hamiltonian K are given by 



I.q 



where H is written as a function of / and r only. 
Now we define J by 



r-2jr 

2tt Jo 



J(K,t) = ±- I Idd> (30) 



where the integral is to be taken round the curve (assumed closed) in the (J, </>) plane given by K = 
constant. If we write I as a function of H and r given by I(H, r), we can use using (|29|l and a Taylor 
series expansion to obtain: 

W,r) -i{k-* {%\, r) = I { K,r) - f (g)^ (§) ^ + 0(f) (31) 
where the final partial derivative is to be evaluated at H = K . Substituting into 130f) . and remembering 



5 More formally we mean that max{|7(t) — /(0)| : t ^ T/e} = O(e) for all one parameter adiabatic variations 
r = r(e, t) with f = 0(e), f = 0(e 2 ) and / evaluated with r = r(e, t). 
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that K is constant in this integral, gives 

J(K,r) - ,(K,r) - *- ^J'' d , + (^) (32) 



u\8t ) I q 2ttu J \3t j I q 

where we have written uj for dH/dl at fixed r. Here we assume that to is bounded away from zero 
throughout the motion. This corresponds to the requirement that the system remains in a quasi-periodic 
state for all times under consideration. Neither to nor (^r) l depend explicitly on the small parameter 
f , so the first result, namely at any given time 

\J-I\=0(f) (35) 

is established. 

Now we have to investigate the time evolution of J. Using that fact that J depends on K and r 
only, we find 

dj /an +(dJ\. = (8j\ Jdj\(dK\ = (bj 



,1/ \0tJ K '\8Kj t K \dt) K + \dK)\dt) IA \dt) ,: ' ( " 

dJ { d\ f r f fdF\ t f 27r fdF\ , . / 

d^ + 0(f 2 ) = 0(t,t 2 ) (37) 



d * \ dt J4,,i\ UJ \ dT Ji. q 2lTUJ Jo \ dT /i, q 

which establishes the second result. 

If the variation in r is adiabatically slow, so that 0(t, f 2 ) = 0(e 2 ), then by integrating dJ/dt from 
to t G [0, T/e] we deduce that \ J(t) — J(0)| = 0(e). Our first result, equation (jSEJ), then implies that 
\I(t) — /(0)| = 0(e), i.e., / is an adiabatic invariant. 



6 Conclusion 

The theory of adiabatic invariants is one of the more confusing aspects of Hamiltonian mechanics in 
undergraduate courses on classical dynamics. An intuitive understanding of the concept is often all that 
is given, and the subtleties ignored. For instance, as we point out, it is not true to say that an adiabatic 
invariant varies much more slowly than the slowly-varying parameter. As we discussed in section|21 the 
action variable typically has a rate of change comparable with that of the slowly- varying parameter. 
However, over suitably defined extended periods of time, the change in the adiabatic quantity is much 
less than the typical change in the quantity being varied. It is in this rather precise sense that the 
quantity is 'invariant'. 

The subtlety in formulating an adequate definition of adiabatic invariance means that a proof of 
the adiabatic invariance of the action variable, for instance, will of necessity require a certain level of 
sophistication. In section [3] we pointed out some of the problems encountered in attempts to prove 
the result directly, i.e., without recourse to the more advanced theory of canonical transformations and 
action-angle variables. Frequently attempts to prove the adiabatic invariance of the action variable 
introduce a time averaging procedure; the average being taken over a time scale over which the system 
is approximately periodic. One may then be tempted to replace the quantity /, which is a geometrical 
construction given by an integral over a region in phase space bounded by a certain constant energy 
contour, with an approximation to I based on the trajectory of the particle in question, a dynamical 
construction requiring the use of Hamilton's equations. This would make no difference if the system 
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were independent of time, but in the time dependent case, one is then left with the problem of keeping 
track of the all the errors made in using these approximations. 

In contrast, our new proof constructs J, a second action variable which has the natural interpretation 
as the phase space average of I over all points on the constant energy contour, or more precisely, over the 
angle variable. This is the analogue of the time-averaging procedure normally used. Perhaps naturally, 
in light of the simplicity with which this fits into the formalism of Hamiltonian mechanics, this leads to 
a far clearer proof of the adiabatic invariance of action variable where error terms can be easily tracked 
at every step of the proof. 

We have used the harmonic oscillator to illustrate our proof and to emphasize the subtleties that 
arise in the theory. For the simplest of variations Hamilton's equations may be solved exactly in terms 
of Airy functions. The asymptotic properties of these functions are well-known and lead to a concrete 
example of the issues we have discussed in this paper. The harmonic oscillator with arbitrary adiabatic 
variation is also simple enough to allow for all of the functions introduced in our general proof to be 
written down in closed form. This provides a useful touchstone for readers following the general proof. 
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